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ABSTRACT 

Gas masses tightly correlate with the virial masses of galaxy clusters, allowing for a precise determination 
of cosmological parameters by means of large-scale X-ray surveys. However, based on recent Suzaku X-ray 
measurements, several groups derived gas mass fractions, / gas , that are considerably larger than the cosmic mean 
at the cluster virial radius, 1?200> questioning the accuracy of the cosmological parameter estimations. Here, we 
use a large suite of cosmological hydrodynamical simulations to study measurement biases of / gas . We employ 
different variants of simulated physics, including radiative gas physics, star formation, and thermal feedback by 
active galactic nuclei, which we show is able to arrest overcooling and to result in constant stellar mass fractions 
for redshifts z < 1 . Computing the mass profiles in 48 angular cones, whose footprints partition the sphere, we 
find anisotropic gas and total mass distributions that imply an angular variance of / gas at the level of 30%. This 
anisotropic distribution originates from the recent formation epoch of clusters and from the strong internal baryon- 
to-dark-matter density bias. In the most extreme cones, / gas can be biased high by a factor of two at R200 in massive 
clusters (M200 ~ 10 15 M©), thereby providing a potential explanation for high / gas measurements by Suzaku. While 
projection lowers this factor, there are other measurement biases that may (partially) compensate. We find that 
at /?200> /gas is biased high by 20% when assuming hydrostatic equilibrium masses, i.e., neglecting the kinetic 
pressure, and by another ~ 10-20% due to the presence of density clumping (depending on mass and dynamical 
state). At larger radii, both measurement biases increase dramatically. While the cluster sample variance of the true 
/ gas decreases to a level of 5% at R200, the sample variance that includes both measurement biases remains fairly 
constant at the level of 10-20% (depending on dynamical state). At the high-mass end (M500 > 2 x 10 14 M o ), the 
true / gas within R500 shows a constant redshift evolution. While this result is in principle encouraging for using gas 
masses to derive cosmological parameters, careful X-ray mocks are needed to control those various measurement 
biases. 

Subject headings: Cosmology: Theory — Galaxies: Clusters: General — Large-Scale Structure of Universe — 
Methods: Numerical 



1. INTRODUCTION 

Galaxy clusters are the rarest and largest gravitationally- 
collapsed objects in the universe, making them very s ensi- 
■ tive t racers of the growth of structure (for a review, see IVoitl 
5_j " 120051) . Clusters exhibit a well-defined number count that 
I steeply falls as mass and redshift increase. The number den- 
sity tail is very sensitive to changes in cosmological parameters 
(as demonstrated by recent m easurements, e.g., by lAllen et al.l 
2008; Vikhlin in et al.l l2009bl) including the dynamical char- 
acteristics of dark energy, primordial non-Gaussianity, and 
the theory of gravity. Hence, future large surveys of clus- 
ters (in the optical, X-ray, and microwave/radio/sub-mm wave- 
lengths) with substantially increased statistics and resolution 
can pote ntially provide a gold mine for fundamental cosmol- 
ogy (e.g.jBattve & Wellerj|2003MHull2003HMaiumdar & Mohri 
2003, 2004; Khede kar et al. 1120 1 OtlGalri et al .11201 21) . However, 
this is only possible if systematics associated with cluster mass 
calibrations and our incomplete knowledge of the physics of 
the intracluster medium (ICM) can be understood and con- 
trolled. As observations have been progressively refined and 
resolutions improved, the clusters have been revealed to be too 
complex for simple sphericalized analytical modelling, calling 
for detailed theoretical work with a necessarily heavy compu- 



tational component. 

The virialized regions of clusters - usually characterized by a 
sphere of radius R 200 enclosing a mean density that is 200 times 
the critical density of the Universe (and higher mean interior 
densities correspond to smaller radii) - are separated from the 
Hubble-flow. However, those regions maintain contact with the 
surrounding filamentary cosmic web through ongoing accre- 
tion and mergers as they evolve. Most of the cluster baryons 
are in the form of a hot (kT ~- (1 — 10) keV) diffuse plasma, the 
ICM, while the remaining baryons are confined to the cluster's 
numerous stars and galaxies. Since cold dark matter (CDM) 
obeys the collisionless Boltzmann equation and since the dy- 
namics of the hot baryons additionally responds to pressure 
forces according to the (collisional) Euler equation (neglecting 
viscosity effects), we expect a strong internal baryon-to-dark- 
matter (DM) density bias. We hope that when averaged over 
cluster scales, the smoothed densities are nearly in the univer- 
sal Hubble-volume-smoothed proportion (or a constant, large 
fraction thereof). But at which radius and cluster scale this 
exactly happens depends on baryonic physics and needs to be 
carefully studied using a combination of hydrodynamical cos- 
mological simulations and ICM observations. 

Using Chandra X-ray observations of a sample of relaxed 
clusters, earlier work has successfully used the gas mass frac- 
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tion in clusters to derive cosmological parameters in agre ement 
with the concordance ACDM model dAllen et al.l 12002, 2004, 
2008). In these analyses, masses were derived assuming spher- 
ical symmetry and hydrostatic equilibrium while allowing for 
a modest constant non-thermal pressure contribution. Similar 
conclusions were obtaine d through a joint X-ra y and Sunyaev- 
Zel'dovich (SZ) analysis (LaRoque et al. 2006) as well as in a 
Chandra sample of X-ray luminou s clusters extending to a red- 
shift of z < 1.3 (Ett oriet al.l l2009). Recently, a more sophisti- 
cated analysis of an X-ray flux-selected sample of 238 clusters 
at z < 0.5, which takes into account various selection effects, 
combi nes gas mass fracti ons, X-ray luminosity, and temper- 
ature dMantz et al.ll2010alah . These authors derive improved 
cosmological parameters that are consistent with ACDM and 
have approximately equally constraining power in comparison 
to other cosmological probes (at the time of publication) on 
the DM density parameter, f2 m , the rms amplitude of the (lin- 
ear) density power spectrum on cluster-mass scales, <7g, and the 
equation of state of dark energy. These data were even used to 
obtain interesting constraints on modified gravity and gen eral 
relativity on cosmological scales (Rapetti et al. 2009, 201CJ) as 
well as on neutrino properties dMantz et al.ll2010bl) . However, 
in order to get unbiased cosmological constraints, the gas mass 
fraction, / gas = M gas /M tot , necessarily needs to be independent 
o f both, redshift a nd mass ( Sa dat et al.l l2005). 

IVikhlinin et alj d2006l) studied gas and total mass profiles 
in a small sample of Chandra clusters of high-quality data 
out to 7?5oo- They find that enclosed / gas profiles within 
/?2500 — 0.4/?5oo are not asymptotically approaching the univer- 
sal baryon fraction according to the cosmic microwave back- 
ground (CMB) observations (with a discrepancy of a factor of 
1.5-2). This discrepancy is becoming smaller for increasing 
cluster masses and radii, but is still noticeable at /?5oo enforc- 
ing the need to understand the behavior of / gas at even larger 
radii. 

Recently the launch of Suzaku enabled to measure the clus- 
ter X-ray emission up to and even beyond /?200 due to its low- 
Earth orbit that ensures a low and stable particle background, 
almost an order of magnitude lower in comparison to previ- 
ous X-ray telescopes. These measurements are critical in im- 
proving our understanding of the formation history of galaxy 
clusters, the dynamical state of the ICM at greater radii with- 
out the need to extrapolate measurements of the core regions of 
clusters, and to gain confidence in assumptions used to derive 
cosmological parameters from cluster surveys. There is now a 
growi ng collection of Suza ku clusters, con sisting of PKS0745 - 
191 dGeorge et alj 120091). Abel l 1795 dBautz et alj 120091) 
Abell 2204 dReiprich et all 120091). Abell 1413 dHoshino et alj 
l20loh. Abell 1689 dKawaharada et al.l l20Toh. Abell 2 142 
dAkamatsu et alJ 120111) . Perse us fsimionescu et al J 120111) . a 
fossil group RX Jl 159+ 5531 dHumphrevet all 120121). Abell 
2029 dWalker et al.ll20ll) . and Hvdra A dSato et alJl2012l) . 

While most of these Suzaku measurements observe an ex- 
cess X-ray emission for radii r > R500 over what is expected 
from cosmological hydrodynamical simulations, there are two 
clusters, Abell 1795 and RX Jl 159+5531, where the emission 
is compatible with the theoretical expectations (within the er- 
ror bars). The excess emission manifests itself in entropy pro- 
files which rise, but less steeply than the predictions of purely 
gravitational hierarchical structure formation (or even start to 
drop beyond 7?2oo)- Alternatively, this implies increasing / gas 
profiles reaching values above the cosmic mean and, in one 
case, even boosting the universal baryon fraction by a factor 
of 1.6 at R 200 when accounting for a stellar mass fraction of 



12% dSimionescu et alJl201 ll) . In some cases, the excess emis- 
sion could be explained by the anisotropic nature of the cos- 
mic filaments connecting to clusters, in particular if there is 
an overlap of the outskirts of two neighboring clusters (as it 
is possibly the case in Abell 2029). Possible physical expla- 
nations include substantial density clumping and/or neglecting 
(respectively underestimating) non-thermal pressure in deriv- 
ing hydrostatic masses. Since the X-ray emissivity is propor- 
tional to the square of the gas density, a clumpy ICM with 
dense regions will induce X-ray surface br ightness fluctuations 
dChurazovet all 120121: ISanders & Fabianl 120121) that enhance 
the emission over the smooth component, biasing the entropy 
lower and / gas higher in comparison to an extrapolation from 
the more homogeneous inner regions. Observational or in- 
strumental reasons for the observed excess emission (which 
is typically a factor of 3-5 below the soft extragalactic X-ray 
background) include fluctuations of unresolved point sources, 
point-spread function leakage from masked point sources, or 
stray light from the bright inner cluster core. 

To better address possible systematics, a large sample of 
clusters with well-defined selection criteria are currently be- 
ing observed, exploiting the strengths of three complemen- 
tary X-ray observatories: Suzaku (low, stable background), 
XMM-Newton (high effective area and sensitivity ), and Chan- 
dra (good spatial resolution) dMiller et al.ll201 ll) . Combining 
Chandra and Suzaku pointings in cluster outskirts increase the 
number of point source detections by a factor of 10 and, hence, 
allow to reduce the uncertainty on the X-ray surface brightness 
by a factor of ~ 3. Early results from this program on Abell 
3378 show entropy and / gas profiles that are in fact compati- 
ble with theoretical expectations. This result is confirmed by 
a stacking analysis of ROSAT clusters t hat detects a steepe ning 
of the density profiles beyond ~ R5Q0 dEckert et alj|201 ll) . in 
agreement with simulation results that include ra diative physics 
and/or account for gas clumping (Nag ai & Laull201 H) - appar- 
ently contradicting recent Suzaku results. 

Hydrodynamical simulations of the formation of clusters 
that include radiative cooling, star formation and stellar feed- 
back predict stellar mass fractions in excess of what is ob- 
served, a decreasing gas mass fraction with red shift, and over- 
predi c t the X-ray luminosity on group scales ( Kravtsov et alj 
2005; Ett ori et al.l 120061 iRoncarelli et all I2006I) . This prob- 
lem is substantially mitigated when includ ing energetic feed- 
back b y active gal actic nuc lei (A G N) (Siiacki & Springel 
1 20061: ISiiacki et alJ 120071 120081: iPuchwein et alj 12 008: 
Booth &Schave 2009: iFabian et al.l I20T0I: iMcCarthv et al] 
20 lOt iDubois et al.l 12010: Tevs sier et all 1201 It I Young et alj 
20111) . Interestingly, including AGN feedback brings the mod- 
eled SZ power spectrum into agreement with data while purely 
radiative models overproduce the measured power on scales 
aroun d 3 arcmin corresponding to angular multipoles of £ ~ 
3000 dBattaglia et alJ20iq|2^11btlShaw"et al.l2010tlTrac et alJ 
1201 ll) . Most notably, the SZ power spectrum of models that 
only account for shock heating and neglect radiative physics 
and hence the condensation of baryo ns into stars produce s too 
much power on scales of I ~ 3000 (Battaglia et al.ll2010h . in 
strong tension with the measurements by ACT and SPT While 
the gas mass fractions in these models are higher than in the 
corresponding radiative models (with and without AGN feed- 
back), they are somewhat lower than the universal baryon frac- 
tions. This demonstrates that for consistency reasons with the 
SZ power spectrum, the Suzaku-inferred high values of / gas 
along some radial arms either are not representative for the en- 
tire virial regions of those clusters or appear to be biased high. 
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Hence this emphasizes that a successful model of the cluster 
physics not only needs to match stellar and gas mass fractions 
as a function of radius and cluster mass, but also has to adhere 
to constraints from SZ surveys, namely scaling relations and 
power spectrum measurements. 

This is the third paper of a series of papers addressing the 
cluster physics of SZ and X-ray surveys. In the first two 
companion papers, we thoroughly scrutinized the influence of 
feedback, non-th ermal pressure and c luster shapes on Y—M 
scaling relations (Batta glia et al.ll201 lal hereafter BBPS 1) and 
the thermal SZ power spectrum (Battaglia et al.ll201 lbl here- 
after BBPS2). The forth paper details the physics of density 
and pressure clumping due to infalling substructure s, accom- 
panying the growth of clusters (Battaglia et al. 2012a, hereafter 
BBPS4), and the fifth provides an information theoretic view 
of clus ters and their non-equilibrium entropies (Battaglia et al.l 
I2012U hereafter BBPS5). 

In this work, we are motivated by current and next gener- 
ation surveys in the X-rays (eROSITA) and those employing 
the Sunyaev-Zel'dovich effect (Atacama Cosmology Telescope 
(ACT), South Pole Telescope (SPT), Planck) to critically re- 
assess the impact of structure formation and physics on / gas 
measurements, in particular in the outskirts of galaxy clus- 
ters. Using a large suite of cosmological simulations, we first 
quantify the bias of / gas measurements on account of non- 
equilibrium processes in the cluster outskirts while allowing 
for cluster sample variance. Then, we provide an extensive 
study of / gas in our various physical models to obtain a solid 
theoretical foundation for cluster cosmology using gas mass 
fractions. In Section|2] we introduce our simulations and mod- 
eled physics. In Section [3] we study how typical assumptions 
in measuring gas fractions, e.g., adopting hydrostatic mass esti- 
mates and neglecting density clumping, bias the measured val- 
ues. In Section|4] we study how / gas varies with radius, cluster 
mass, redshift, and simulated physics, and compare our simu- 
lation values to observational data. We conclude in Section[5] 

2. COSMOLOGICAL SIMULATIONS AND CLUSTER DATA SET 

We use smoothed particle hydrodynamic (SPH) simulations 
of large-scale periodic boxes which provide us with large clus- 
ter samples to accurately characterize ICM properties over 
large ranges of cluster masses a nd redshifts. W e use a modi- 
fied version of the GADGET- 2 dSpringelll2005l) code. Our se- 
quence of periodic boxes have lengths of 165 and 330 /T 1 Mpc 
filled with A^dm = Ngas = 256 3 and 512 3 dark matter (DM) and 
gas particles, respectively. This choice maintains the same 
initial gas particle mass m gas = 3.2 x IO^/t'M©, DM particle 
mass »jdm = 1-54 x 10 10 h~ l M@, and a minimum gravitational 
smoothing length e s = 20h~ l kpc; our SPH densities were com- 
puted with 32 neighbours. For our standard calculations, we 
adopt a tilted ACDM cosmology, with total matter density (in 
units of the critical) Q m = JIdm + = 0.25, baryon density fit, 
= 0.043, cosmological constant Q\ = 0.75, Hubble parameter 
h = 0.72 in units of 100 km s -1 Mpc -1 , spectral index of the 
primordial power-spectrum n s = 0.96 and a% = 0.8. 

We show results from simulations employing three differ- 
ent variants of the the simulated physics: the first only ac- 
counts for gravitational shock heating; the second additionally 
accounts for radiative cooling, star formation, supernova feed- 
back, galactic winds, and cosmic rays physics and is referred to 
as radiative cooling model; and the last accounts for radiative 
physics (radiative cooling, star formation, supernova feedback) 



and thermal feedback by AGNs and is referred to as AGN feed- 
back model. Radiative cooling and heating were computed as- 
suming an optically thin gas of a pure hydrogen and helium pri- 
mordial composition in a time-dependent, spatially uniform ul- 
traviolet background. Star formation and supernovae feedback 
were modelled usin g the hybrid multiphase mod el for the inter- 
stellar medium of Springel & Hernquist (2003). The CR pop- 
ulation is modelled as a relativistic population of protons de- 
scribed by an isotropic power-law distribution f unction in mo- 
mentum space with a spectral index of a = 2.3 {Pfrommer et alj 
120061 !EnBlin1it1d1l2007t Uubelgas et alJl2008h . With those pa- 
rameters, the CR pressure modifies the SZ effect at most at the 
percent level and cause s a reduction of the r esulting integrated 
Compton-y parameter ( Pfrommer et al. 2007). The AGN feed- 
bac k prescription includ ed in the simulations (for more details 
see lBattaglia et al.l l2010) allows for lower resolution and hence 
can be applied to large-scale structure simulations. It couples 
the black hole accretion rate to the g lobal star format i on rat e 
(SFR) of the cluster, as suggested by Thompso n et al.l ((2005). 
The thermal energy is injected into the ICM such that it is pro- 
portional to the star formation within a given central spherical 
region. 

We use a two-step algorithm to compute the virial mass of 
a cluster. First, we find all clusters in a given snapshot using 
a friends-of-friends (FOF) algorithm. Then, using a spherical 
overdensity method with the FOF values as starting estimates, 
we iteratively calculate the center of mass, the virial radius, Ra, 
and mass, Ma, contained within Ra, and compute the radially 
averaged profiles of a given quantity with radii scaled by Ra- 
We define Ra as the radius at which the mean interior density 
equals A times the critical density, e.g., for A = 200 or 500. 
Our cluster sample (for each simulation type) consists of ~ 
1300 clusters with M 2 oo > 7 x 10 13 M and ~ 800 clusters with 
M 2 oo > 1 x 1O 14 M atz = 0. 



3. BIASES OF GAS MASS FRACTION MEASUREMENTS 

We are interested in the enclosed gas mass fractions that we 
define as the ratio of gas mass and total mass inside a given 
radius, i.e., / gas (< r) = M gas (< r)/M lot (< r). All our values 
of / gas and the stellar mass fraction, / star , are normalized by 
the cosmic baryon fraction / c = Q&/ ^ m that should be substi- 
tuted by the current best fit value from cosmological precision 
experiments and was assumed to be f c = 0.172 in our simu- 
lations, which compares well to the currently favored value 
of 0.166 (Komatsu et al] |20Tl1) . Here, we study how typical 
assumptions in measuring gas fractions, e.g., adopting hydro- 
static mass estimates and neglecting density clumping, bias the 
measured values. We additionally address how the sample vari- 
ance of clusters and the angular orientation variance within 
clusters impacts inferred values for / gas , i.e., how an / gas mea- 
surement from a mosaic X-ray observation along a radial arm 
varies with angular orientation depending on whether that arm 
intersects a cosmic filament or not. 

In Figure[U we show the median enclosed gas mass fractions, 
/gas(< r) as a function of radius for different cluster masses in 
our AGN feedback model; on the left panel for our sample of 
all clusters, and on the right panel for the 3rd most relaxed 
clusters. To define a relaxed cluster, we split the population 
of clusters into lower, middle, and upper bands of the ratio of 
kinetic-to-thermal energy, K/U, within R200, which we use as 
a measure of dynamical state for the clusters. We define the 
internal kinetic energy, K, and thermal energy, U, of a cluster 




Figure 1. The assumption of hydrostatic equilibrium and of a smoothly distributed ICM introduces a bias in the measurement of / gas . We show the median radial 
profiles of the gas mass fraction, /gas(< r) = M gas (< r)/Mt t(< r), normalized by the universal baryon fraction for different assumptions in calculating M gas and 
M tot as a function of cluster mass. We compare the true gas mass fraction, / gas .true (solid), a biased X-ray measurement of / gas due to the assumption of HSE, i.e., 
without accounting for kinetic pressure support, /„ a s.HSE (dotted), as well as a biased measurement of / gas due to the assumption of HSE and additionally neglecting 
density clumping, / gils .HSE+clump (dashed). The profiles in the left panel reflect the median of the distribution of all clusters and the right panel shows the median of a 
sub-sample of the third most relaxed clusters (with the lowest ratio of total kinetic-to-thermal energy, K/U). There is a bias in / gas introduced by the assumption of 
hydrostatic equilibrium even for relaxed clusters. Here, the pressure derivative was smoothed with adjacent radial bins. 



as 

K«R 20 o)^J2 3mS 7 Pklai ' W 

i 

U«R 200 )^ 3 -^^l, (2) 

where m and p are the gas mass and the SPH density, respec- 
tively for all particles i less than /?200, Pth = kTp/(fj,m p ) denotes 
the thermal pressure (where T and fi denote the temperature 
and mean molecular weight, respectively), and the kinetic pres- 
sure is defined as 

P kin = pSv 2 /3, (3) 
Sv = a (v-v) + aH(a)(x-x) . (4) 

The code employs comoving peculiar velocities, which are 
translated into the internal cluster velocities relative to the over- 
all mean cluster velocity in the Hubble flow by the relation 
given, where H(a) is the Hubble function, a is the scale factor, 
v (= dr / df ) is the peculiar velocity, and x is the comoving po- 
sition of each particle. The gas-particle-averaged cluster bulk 
flow within P200 is v and the center of mass within /?200 is x. 
Since the most massive clusters are on average also the most 
disturbed systems, only two out of 25 clusters in our most mas- 
sive cluster bin (M200 > 5.1 x 10 14 Mq) meet our criterion of 
being a relaxed cluster, i.e., are among the lowest third when 
selected by the ratio of kinetic-to-thermal energy, K/U. Hence, 
in order to improve the statistics of our most massive bin of re- 
laxed clusters, we decided to do a 3-split on K/U in the largest 
cluster mass bin only. 



As shown in Figure [T] / gas increases with radius and ap- 
proaches values within 2i?200 that range from 0.8 to 0.85 (in 
units of f c ) for the mass interval from M200 — 1O 14 M to 
5 x 10 14 Mq. Smaller systems show smaller values of / gas at 
every radius owing to the shallower potential wells in compar- 
ison to bigger systems. In groups and small clusters, AGN 
feedback is able to push the gas well beyond R200 or halt its 
initial infall. In the process of hierarchical structure formation, 
these gas-poor groups merge to form larger clusters. Owing to 
their deeper potentials, those larger clusters provide a greater 
gravitational pull on the gas that is sitting beyond 7?2oo- As a 
result, / gas increases at e ach radius in larger clusters. (We will 
show later on in Section l4~3l that the underlying assumption of 
a constant / gas as a function of redshift at fixed cluster mass is 
justified in our simulations.) 

These simulated values of fg as /f c < 0.8 for all cluster masses 
are substantially smaller than Suzaku estimates of / gas //r ~ 1 .4 
or even 1.6 within P200 if one accoun ts for a stellar mass frac- 
tion of 12% (Simi onescu et alj|201 lh . Two sources that could 
potentially bias the measurement of gas fractions have been 
suggested, the assumption of hydrostatic masses and neglect- 
ing density clumping due to infalling (sub-)structures follow- 
ing the ongoing growth of clusters. Hence we will discuss each 
of these effects in turn. 

3.1. Hydrostatic mass bias 

Even if clusters are in hydrostatic equilibrium, balancing the 
gravitational force to the pressure gradient yields 

VP = pg — > -pGM(< r)r/ r 2 for spherical symmetry. (5) 
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Assuming that all the pressure in Eq. (O is thermal (P = f t h) 
gives an incorrect estimate for the mass since the non-thermal 
press ure, in particular the kinetic pressure has to be included 
(e.g.. lEvrarcMf990tlRasia et alj2004lLaul^l2009l) . For clar- 
ity we define Mhse to be the mass derived using P = P±. 

Comparing M H se to the true mass inside a given radius, 
Mhse on average underestimates M to t by 20-25% within Z?200 
(BBPS1) and biases / gas accordingly high (dotted lines in Fig- 
ure [T). This Mhse bias is almost independent of cluster mass 
out to R200 as can be seen by the fractional difference of / gas 
(lower panels of Figure [TJ. However, individual clusters can 
stray from this average deviation, since each cluster has a 
unique dynamical state and formation history. This is sug- 
gested by a scatter of ~ 5% betwe en the 2 5 th and 75 th per- 
centiles of the complete distribution (BBPS1). For approaches 
that assume a constant Mhse bias, the fractional difference of 
/gas.HSE would be accordingly shifted downwards (by the same 
factor in the lower panels of Figure [T). 

When we select a sample of the third most relaxed clus- 
ters, the correction factor will necessarily be smaller. We find 
that the bias on / gas for this sample is reduced to 15-20% 
within /?200- In a small sample of zoomed-cluster simulations 
that has been selected against major mergers the hydrostatic 
mass correction was foun d to be of the order M H se ~ 10- 
15% (iKravtsov et alJl2006h . However, in both cases, we see 
that outside R200, the hydrostatic bias rises steeply as i\in/Ah 
approaches large values and eventually reaches equipartition 
around 2R 200 (at z = 0) with the exact equipartition radius de- 
pending on cluster mass. 

3.2. Clumping bias 

The emissivity of thermal bremsstrahlung scales with the 
square of the gas density. Hence the X-ray-inferred gas density 
profile is biased high by the presence of a clumped medium and 
attains a bias of the square root of the density clumping factor, 
that we define by 



Averages are taken to be volume averages. We also define the 
SPH volume bias, 

i i 

that we expect to be identical to unity everywhere. While this 
is the case for small radii, the SPH volume bias deviates pro- 
gressively for larger radii. It reaches a negative bias of 30% at 
3^200 (z = 0), almost independent on the simulated physics and 
cluster masses, but not of redshi ft since it doubles at z = 1 when 
radii are scaled to 1?200 (BBPS4). To obtain unbiased estimates, 
we decided to redefine the clumping factor, yielding 

Ci.p = C2. P Co, (8) 

where the multiplication with the volume bias factor, Co, has 
been chosen to cancel the SPH volume in the weighting. Inside 
R200, the profile of \JC2.p increases moderately with radius but 
remains < 1 .2 for r < R200 f° r all clusters. Outside R20Q, it starts 
to rise sharply and reaches values of 10 at r ~ 3/?200 and 4i?200 
for large clusters and groups, respectively, almost independent 
of the simulated physics (BBPS4). At each radius, there is a 
clear trend of an increasing clumping factor for larger clusters, 
which have on average a higher mass accretion rate, i.e., they 



are assembling at the current epoch and hence are dynamically 
younger. 

Cold, dense gas clumps that are ram-pressure stripped from 
galaxies and find themselves in pressure equilibrium with the 
ambient hot ICM add an artificial bias to the density clumping 
factor. This bias needs to be accounted for when comparing 
to X-ray observations by employing, e.g., a temperature cut. 
Hence, in our calculations of the density clumping factor, we 
only consider SPH particles with T > 10 6 K that can be ob- 
served by their X-ray emission. 

To estimate the clumping bias of the enclosed gas mass frac- 
tions, we compute a clumping-biased density profile for each 
cluster by means of PciumpW = vZ&.pW Arae( r )> an d integrate 
this to obtain a clumping-biased gas mass profile, 

M gas , dump (< r) = 4-n / drr 2 p clump (r). (9) 
Jo 

Here, we only asses the effect of clumping bias on M gas . In 
principle, density clumping should also impact the determina- 
tion of M tot . However, there are different subtle effects that 
may partially compensate. Those include larger Poisson scat- 
ter in density profiles implied by clumping, which then have a 
noisier derivative, calling for some smoothing procedure which 
may bias the inferred masses. Also, the X-ray emissivity is en- 
hanced by the presence of self -bound (dense) structures which 
increases the logarithmic derivative of the density and hence 
M tot . Conversely, these dense substructures have a lower tem- 
perature. Fitting a single temperature spectrum to the data will 
then bias the inferred temperature and hence M tot low. Simu- 
lated X-ray observations are needed to carefully address these 
points that shall be subject to future work. 

Clumping bias increases inferred values for / gas - similar to 
the Mhse bias. In Figure Q] we show the effect of both bi- 
ases that we denote by / gas .HSE.ciump (dashed lines). As a re- 
sult of both biases, / gas realizes its universal value at consider- 
ably smaller radii in comparison to the true / gas . In our sample 
of relaxed clusters, the universal value is reached outside /?2oo 
whereas in our all-cluster sample, these high values are already 
obtained at 0. 57? 200 for big clusters with M200 > 5 x 10 14 M Q . 
Interestingly, the biased values for / gas increase sharply outside 
the virial radius, mainly driven by the clumping bias for large 
systems and with a substantial contribution of Mhse bias for 
smaller clusters. / gas is biased high by a factor of two at 27?2oo 
for large clusters and at 3/?200 for groups (see FigureQ]i. 

What is the relative contribution of these measurement bi- 
ases as a function of radius and cluster mass? In our all-cluster 
sample, clumping bias has a small effect at the cluster center 
but is continuously increasing towards Z?20Q. At R200, it ap- 
proximately contributes half as much as the Mhse bias to the 
overall bias with an absolute bias of 10-20%, depending on 
cluster mass. In our relaxed cluster sample, the clumping bias 
remains even more subdominant in comparison to the Mhse 
bias. However, in both cluster samples, the clumping bias in- 
creases dramatically beyond /?20o with an increasing bias with 
cluster mass. 

3.3. Cluster sample variance 

It has been suggested that sample variance of clusters could 
possibly account for values of / gas larger than the cosmic mean, 
i.e., when a cluster experiences an ongoing merger or has close- 
by cluster companions. To address such questions, our large 
sample of clusters is ideal. It consists of 800 clusters with 
M 2 oo > 1 x 1O 14 M and has 200 clusters with kT spec > 2.5 keV 
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All clusters: Relaxed clusters: 




Figure 2. Cluster sample variance of /gas. Upper panels: we show the gas mass fraction / gas for different mass intervals (lines indicate median values, shaded regions 
are bounded by percentiles such that they contain 68.3% of the probability, centered upon the median). Lower panels: we show the relative cluster-to-cluster scatter 
of the gas mass fraction, oy , that we define by the difference of the upper and lower (1 — it) percentiles normalized by twice the median of / gas . Here, we compare 
radial profiles of the true gas mass fraction, / gas ,true (solid), a biased X-ray measurement of / gas due to the assumption of HSE, i.e., without accounting for kinetic 
pressure support, / eas hse (dotted), as well as a biased measurement of / gas due to the assumption of HSE and additionally neglecting density clumping, / gas .HSE+clump 
(dashed). The profiles in the left panel include all clusters and the right panel shows the median of a sub-sample of the third most relaxed clusters (with the lowest 
ratio of total kinetic-to-thermal energy, K/U). While the scatter of / gas ,tme decreases for larger radii, the scatter of the biased / gas ,HSE+clump remains fairly constant at 
a level of oy //gas — 0. 1 —0.2 (depending on the dynamical state and cluster mass) within K200 and even stalls to increase for larger radii due to the strong clumping 
term. (We do not show the shaded regions of / gas .HSE for visual clarity, but show this information in the lower panels.) 
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Figure 3. Upper panels: we show the profiles of A/ gas (left), Mtot (middle), and / gas (right) in 48 angular HEALPix cones (directed such that they subtend an equal 
area pixelization on the sphere), centered on one cluster of mass M200 ~ U2 X 1O 15 M0. (Note that we only show every other profile for visual purposes.) We 
overplot the median (solid black) as well as the upper and lower (1 — a) percentiles (dotted black). Lower panels: we show the variations of M gas , Mtot, and / gas in 
the different cones relative to their corresponding median. Note the different scalings at the y-axis (logarithmic for the mass profiles, linear scale for / gas ). There is 
considerable scatter as well as the presence of a few enormous outliers in the M gas and M to t cones - depending on the angular direction - which is reduced for / gas . 
This implies spatial correlations of the DM and gas distribution through self-bound substructures that are still holding on to (some) of their gas. 



(corresponding to M 50 o > 2.15 x 10 14 M Q ). In Figure |2] we 
compare the cluster-to-cluster scatter of / gas in a sample of all 
clusters and a sub-sample of the third most relaxed clusters 
(as defined by the third of all clusters with the lowest ratio of 
kinetic-to-thermal energy, K/U, within R200, see Equations ([TJ 
and (0). The sample variance of the true / gas decreases with 
increasing radius from a relative 1-er scatter of 20% (10-15%) 
at 0. l/?200 for our sample of all (relaxed) clusters down to ~ 5% 
at 1?200- This demonstrates the smaller influence of AGN feed- 
back on large scales especially for big clusters where the grav- 
ity of the DM potential starts to dominate. In contrast, the 



scatter of the biased quantity / gas ,HSE remains fairly constant 
at a level of o/„ as // gas ~ 0.1 -0.2 (depending on the dynami- 
cal state and mass of clusters). Within R200, the sample vari- 
ance / gas due to the clumping bias is negligible, but it starts to 
rise steeply outside /?200> the innermost locations of accretion 
shocks (along filaments) that dissipate kinetic accretion energy 
and smooth the density distribution interior to this radius. 

The 1 - 0- envelope of / gas ,HSE+ciump starts to rise dramatically 
for r > 7?2oo for our all-cluster sample. This rise is shifted to 
somewhat larger radii for our sample of relaxed clusters (Fig- 
ure O. However, at R 200 we obtain an upper 1 - it percentile 
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Figure 4. Top panels: we compute the variance of the M gas profiles across 48 angular HEALPix cones (dashed lines in Figure[5J for each cluster and show the mean 
of this variance across all clusters (left) and across a sub-sample of the third most relaxed clusters with the lowest ratio of kinetic-to-thermal energy, K/U (right). We 
define the angular variance of A/ gas , &M ai!i , by the difference of the upper and lower (1 — cr) percentiles normalized by twice the median of M gas . We compare the mass 
measurement (dashed) and our estimate of the intrinsic angular variance (solid) where a Poisson noise term due to the finite particle number within the cones has 
been subtracted in quadrature. Note that this over-corrects the intrinsic variance at small radii as the Poisson noise estimate caps the noise contribution owing to the 
true sub-Poissonian distribution of SPH particles, which get rearranged as a result of the mutual repulsive hydrodynamical forces. Middle and bottom panels: same 
as above, but for M tot and / g!is , respectively, where M tot is calculated by summing up the mass distribution of the gaseous, stellar, and DM components in those cones. 
Different colors correspond to different cluster mass ranges (indicated in the legends). Due to its collisional nature, gas is able to dissipate part of its kinetic accretion 
energy more efficiently than DM, thereby reducing its angular variance, in particular for relaxed systems (on the right). This is in contrast to the collisionless DM 
component (dominating the total mass) which shows a larger angular variance of its mass distribution. 
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of /gas,HSE+cium P for our massive clusters (M200 > 3 x 1O 14 M ) 
that is 30% above the universal value for our all-cluster sam- 
ple and only 10% for our subsample of relaxed clusters. Clus- 
ter variance alone remains unlikely to support values that are 
1.6 times the cosmic mean - the measured value at R2Q0 in 
Perseus when accountin g for a stellar mass fraction of 12% 
(Simio nescu et al.11201 U even when accounting for the mod- 
erate increase of this number by ~ 10% for the most massive 
clusters around M200 ~ 1O 15 M ). Moreover, clusters that are 
lying outside the 1 — a interval are either merging systems or are 
caught in a pre-merger state, properties that can easily be veri- 
fied through optical or X-ray observations of the local environ- 
ment of such a cluster (unless the merger axis is aligned with 
the line-of-sight, but such a configuration is rare and should be 
easily separable in a statistical sample of clusters). 

3.4. Angular variance of the mass distribution within a single 

cluster 

In nearby clusters, the large angular extent of clusters ren- 
ders it impossible to cover the entire virial region of a clus- 
ter with X-ray pointings. Instead, a limited number of mosaic 
pointings along radial arms are being observed which opens 
the possibility of obtaining a biased answer for / gas if there 
is a considerable angular variance. Two effects may play a 
role. First, a large mass fraction is being accreted along fila- 
ments. These anisotropic ICM features associated with large- 
scale structure environment suggest that the thermalization of 
the ICM occurs at smaller radii along overdense filamentary 
structures (that exert a greater ram pressure on the hot ICM) 
than along low-density void regions. Second, since there is 
a strong internal baryon-to-dark-matter density bias, a con- 
sequence of collisional-to-collisionless physics, this could in 
principle also lead to a significant angular variance depending 
on the differing efficiencies of gas thermalization and the re- 
laxation processes of the DM component. 

To quantify these considerations, in Figure [3] we show the 
profile of Mgas, M lot , and / gas in angular cones centered on one 
massive cluster of M200 w 1.2 x 10 15 M Q . We chose the an- 
gular directions of these 48 cones such that they subtend an 
equal area pixelization on the sphere (follow ing the HEALPix 
scheme with N & \d e = 2; Gorski et al.l [2005). The number of 
cones is a compromise between obtaining high-enough angu- 
lar resolution to break the angular smoothing that is necessarily 
present in sphericalized profiles and to resolve individual sub- 
structures and cosmic filaments in these cones. At the same 
time, we aim at minimizing scatter due to our finite number of 
particles. To calculate M tot , we sum up the mass distribution of 
the gaseous, stellar, and DM components in those conesQ 

We observe considerable scatter in all three distributions, 
M gas , M tot , and / gas , depending on the angular direction. There 
are a few enormous outliers in the M gas and M tot cones, deviat- 
ing upwards up to a factor of 3 from the median (shown with 
solid black) in the case of M„ as and even beyond for M tot . The 
right-most panel in Figure [3] shows that in the most extreme 
cases, the cone profiles of / gas can be biased high by a factor 
of two in such a massive cluster with M200 ~ 1O 15 M , poten- 
tially providing an explanation of the boos ted values of f gas 
repor ted in the north-western arm in Perseus (Simionescu et al. 
1201 lh . Care should be taken when interpreting these / gas val- 

1 For the gas component, we sort the volume fraction of an SPH smooth- 
ing kernel into the corresponding radial bins, but neglect the extension of 
the smoothing kernel in angular directions for performance reasons (see Ap- 
pendix [a] 



ues in cones, since projection effects may lower these values. 
However, as demonstrated in Figure[T] hydrostatic mass bias as 
well as density clumping may (partially) compensate for these 
projection effects. 

Only in the M lot cones, those outliers can bee seen down 
to small radii r < O.I-R200 which is expected as DM substruc- 
ture only obeys collisionless physics that allows substructures 
to penetrate deep into the halo without too much affecting the 
dense central subhalo structures by tidal stripping. In contrast, 
the gaseous content has been removed by means of shocks or 
ram-pressure stripping at small radii. These large outliers are 
absent in the distribution of / gas cones which implies a smaller 
variance (dashed black). This suggests that these outliers are 
caused by self -bound substructures that are still holding on to 
(some) of their gas at larger cluster radii, i.e., there exists a tight 
spatial correlation of the DM and gas distribution. 

3.5. Variation of the mass distribution with angle across 
cluster samples 

In order to obtain statistically valid answers on how this scat- 
ter across different cones varies among our cluster sample, we 
first compute the variance of the profiles of M gas , M tot , and 
/ gas across 48 angular cones (dashed lines in Figure [3]) for all 
clusters. Since the distribution within a single cluster is non- 
Gaussian and positively skewed (by the presence of substruc- 
tures), we define the angular variance of M gas , er^ , by the dif- 
ference of the upper and lower 1 — a percentiles, normalized by 
twice the median of M gas (shown in the top panels of Figure [4] 
and color coded by cluster mass). The same quantity is being 
calculated for M tot and/ gas (lower panels of FigureE} across all 
clusters (left) and across a sub-sample of the third most relaxed 
clusters (right). 

As shown in Figure [4] there is an upturn of the raw vari- 
ances (dashed) visible at small radii and in particular for galaxy 
groups. This is caused by the noisy distribution of a compara- 
bly small number of SPH particles in the central cluster regions 
in combination with the smaller radial bin size of our logarith- 
mic equally-spaced grid towards the center, when distributing 
the particle masses among the angular cones. As a result of the 
mutual repulsive hydrodynamical forces of neighboring SPH 
particles, they maintain a sub-Poissonian distribution with re- 
pulsive short distance correlations on scales of order the cor- 
responding inter-particle separations. The full characterization 
of such a distribution with dynamical evolution of SPH is not 
possible in terms of known distributions. Hence we decided 
to cap this numerical noise contribution by treating it as Pois- 
sonian and subtract a Poissonian noise term in quadrature from 
the raw angular variance. We caution that this over-corrects the 
intrinsic variance at small radii, causing an artificial decrease 
of the bin-noise subtracted variance (solid). We emphasize that 
our focus here is on large scales with r > (0.2— 0.5)/?2oo (where 
the range corresponds to large and small clusters, respectively) 
are not subject to this numerical limitation. Hence in order to 
improve upon our analysis in the central regions, an increase in 
numerical resolution would be needed. 

The angular variance of M gas does not strongly depend on 
mass and reaches values around 30% at /?200 with a moderate 
increase at larger radii. In contrast, the angular variance of M tot 
is steadily increasing and reaches values around 45 - 50%, de- 
pending on cluster mass (with the upper range indicating larger 
clusters that are on average dynamically younger). These val- 
ues are only moderately reduced by ~ 5 percentage points for 
our relaxed cluster sample. The different angular variances of 
M gas and M tot can be understood by the collisional nature of the 
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gas that is able to dissipate part of its kinetic accretion energy 
more efficiently (through shocks and ram-pressure stripping of 
the gaseous contents of substructures), thereby reducing its an- 
gular variance. This is in contrast to the collisionless DM com- 
ponent (dominating the total mass) which shows a larger angu- 
lar variance of its mass distribution. 

Interestingly, the Poisson noise-corrected angular variance 
of / gas remains almost constant with radius at the level of 30- 
35% not only within the clusters but also for r > R200 where 
the angular variances of the gas and total mass are increasing. 
Furthermore, ay // gas depends very weakly on cluster mass. 
Both of these properties can be understood by the fact that at 
these large radii, most of the accreting substructures have not 
yet passed the accretion shock nor are they affected by ram- 
pressure stripping as they still move through the dilute warm- 
hot inter-galactic medium. Hence these self-bound substruc- 
tures are still holding on to their gas, implying large spatial 
correlations of the DM and gas distribution (see also Figure[3]l. 

The most striking feature of all our Poisson noise-corrected 
angular variances is the rather constant behavior with radius, 
in particular for the angular variances of M gas and / gas , as well 
as the moderate radial increase of that of M tot . In principle, the 
mass distribution in cones allows for finer sampling of the clus- 
ter shape beyond the moment-of-inertia tensor and sensitive to 
higher-order moments - a measure that would be similar to 
density clumping for sufficiently high multipoles. In practice, 
this is limited by the Poisson noise due to our finite particle res- 
olution. Moreover in Figure [4] we concentrated on the radial 
profiles of the second moment of the angular variance which 
is intimately related but not identical to cluster ellipticities as 
obtained through eigen values o f the moment-of-inertia tensor 
as a function of radius (BBPS1). The ellipticity of the gas den- 
sity also shows a constant behavior inside 7?2oo while it is rising 
at larger radii. In contrast, the DM ellipticity shows already a 
moderate increase with radiu s which only starts to rise more 
steeply outside /?200 (BBPS1) - very similar to our findings for 
the angular variance of M tot that is dominated by DM. 

4. GAS MASS FRACTIONS IN COSMOLOGICAL SIMULATIONS 

Here, we present results on how / gas varies with simulated 
physics, radius, cluster mass, and redshift. 

4.1. How different physics determines / gas and / star 

There is considerable uncertainty about the correct physics 
governing the evolution of the ICM. Hence it is instructive to 
compare how differently simulated physics impact the radial 
profiles of / gas (Figure [5]) and / star (Figure [6j. This is supple- 
mented by a figure demonstrating the cluster mass dependence 
of / gas and / stal within R^q and /?200 for differently simulated 
physics (Figure|7). 

In our shock heating-only model, / gas approaches a value of 
0.93 ± 0.03 in units of f c at r ~ Z?200 (which is reduced by one 
percentage point for larger clusters). Outside the virial radius, 
/ gas is continuously rising to asymptote to the cosmic mean 
only at 4R 200- This may be partially due to a transfer of energy 
from DM to the gas during group/cluster mergers. In an equal 
mass collision, SPH simulations show that the DM can transfer 
up to 10% of its energy, and even a 10:1 merge r may be able to 
reduc e the DM energy by a couple of percent dMcCarthy et alJ 
120071) . There is also the possibility that some of this gas defi- 
ciency may be related to the SPH technique that we use here. 
There appears to be a discrepancy to adaptive mesh-refinement 
simulations that find higher values of /& ~ 0.97 ± 0.03 within 
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Figure 5. Different levels of sophistication of the simulated physics of the 
ICM imply different profiles for / gas (top panel). We show / gas (< r) normal- 
ized by the universal baryon fraction for different physics models at two red- 
shifts, z = (solid) and z = 1 (dash-dotted). At z = 0, we also show / gas (< r) 
assuming hydrostatic equilibrium (dashed). In the bottom panel, we show the 
differences with respect to the true / gas at z = 0. There is very little redshift 
evolution between z = and z = 1, in particular for large radii and for our AGN 
feedback model. The bias in /g as when assuming hydrostatic equilibrium is 
almost independent of the simulated physics. 
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Figure 6. Whether or not we include AGN feedback does affect the me- 
dian profile of the stellar mass fraction, /star- We show / s tai(< r) = M star (< 
r)/Mtot(< r) normalized by the universal baryon fraction at z = (solid) 
and z. = 1 (dash-dotted). While ongoing cooling increases the central stellar 
mass fraction to unphysical high values in our model with cooling and star 
formation-only, AGN feedback is able to self-regulate the core region and to 
suppress the forming cooling flow. 



/?vi r (Kravtsov et al. 2005, see also Section |431 for a detailed 
discussion). At small radii, the shock-heated ICM provides 
pressure that causes the inner slope of the gas density profile to 
be rather flat, w hile the collisionless DM exhibits a cusp (e.g., 
Gao et al]l2012l) . As a result, this causes decreasing gas frac- 
tion inside the DM scaling radius, i.e., for r < 0.2/?2oo- 
Our models with AGN feedback and radiative cooling have 
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Figure 7. We show gas, stellar, and baryon mass fractions, / gas (< r), /star(< r), and /&(< r), (mean and standard deviations), computed within a fixed radius (S500 
and R200X as a function of cluster mass at z = (left panel) and z= 1 (right panel). Different colors indicate different simulated physics models: AGN feedback (red), 
radiative cooling (blue), and shock heating-only (green). Note that f), is almost identical within the scatter for the integration radii R500 and R200- In our cooling and 
star formation model at z = 0, we find an increasing stellar mass fraction towards smaller clusters that is accompanied by a decreasing gas mass fraction, particularly 
within K5oo- This trend of a decreasing / gas towards smaller clusters is amplified in our AGN feedback model, while the stellar mass fraction remains almost constant 
across our mass range probed. AGN feedback pushes gas out of smaller clusters beyond K200 which causes a reduction of the baryon mass fraction in comparison to 
our model without AGN feedback. In contrast, the trends of / slai with cluster mass are almost absent at z = 1 ■ 



similar values of / gas at z = 0, in particular for r > ^500- How- 
ever, in the AGN feedback model, / gas only marginally evolves 
with z, while it increases for higher z in the case of radia- 
tive cooling and star formation. This also manifests itself in 
the constant evolution of the / star profile in the AGN feed- 
back model. In contrast, the stellar mass is increasing with 
time in the radiative model without AGN feedback due to the 
well-known overcooling problem in cosmological simulations 
of galaxy clusters if coolin g is not co unteracted by any feed- 
back process (e.g., Borgani & Kravtsovll2009h . 

For r > R200, the values of / gas in different radiative models 
(as well as for different redshifts) approach each other. How- 
ever, in comparison to our shock heating-only model, there is 
still a gas mass deficit of 10% at 42?200- In the raditive cool- 
ing model, the stellar mass fraction makes up for this deficit 
(Figure |6j while in the case of AGN feedback, stellar mass 
fractions are reduced by a factor of ~ 2 in comparison to our 
radiative cooling model. (Note that / star in our radiative cool- 
ing models depends on numerical resolution as this model is 
physically not converged.) The remaining gas mass has been 
pushed to even larger ra dii through the action of AGN feedback 
(McCarthy et aT] |201 lb . Interestingly, these lower values for 
/gas in our radiative models in comparison to our non-radiative 
models im ply a temperature i ncrease in the outer parts (e.g., 
Figure 4 of iFabian et alj|2010l) that is necessary to balance the 
gravitational cluster potential. 

All our models show approximately a similar relative Mhse 
bias (bottom panel of Figure |5). However, there are some no- 
table differences. At large radii, the radiative cooling model 
exhibits the largest bias, while our AGN feedback model shows 
the lowest bias. This is presumably due to AGN outflows de- 
positing heat at these large radii (early on) which manifests 
itself in an inc reased pressure in co mparison to the radiative 
cooling model (Battagli a et alj|20lol) . This increased thermal 
pressure moves the average location of accretion shocks to 
larger radii (since these form at the location where infalling 
ram pressure balances the thermal pressure of the ICM), which 



causes an deceleration of the accreting material already further 
out, and hence slightly decreases the average kinetic pressure 
support in the AGN feedback model. Hence, the bias correc- 
tions shown in Figure Q] are conservative given that the AGN 
feedback model shows an overall lower bias in comparison to 
the other models. 

4.2. Gas and stellar mass fractions across cluster masses 

How do the gas and stellar mass fractions, / gas and / star , as 
well as the sum of both, the baryon mass fraction, ft,, vary with 
cluster mass? In Figure [7j we show / gas and / star for our 3 
different physics mod els as a function of cluster mass. As dis- 
cussed in Section |4~T1 our model that only accounts for shock 
heating approaches a value of / ga s//c = 0.93 ± 0.03 (1 — er) at 
r ~ R 20l) (which is reduced by one percentage point for larger 
clusters). These lower values of / gas represent the pivot point 
from which our radiative models can reduce / gas furthermore 
through conversion of baryons into stars or by expelling it from 
the cluster region (in the case of our AGN feedback model). 

In our radiative cooling model, / gas increases moderately 
with cluster mass at Rsa) to level off at a value of / gas // c = 0.79 
for big clusters. This trend is weakened at /?200- The oppo- 
site behavior is seen for / star so that the resulting baryon frac- 
tion almost exactly sums to a value that is seen in our shock 
heating-only model. At Z?200> the scatter in / gas drops by 1-2 
percentage points in comparison to the scatter at R500 for both 
of our radiative models. 

While in our AGN feedback model, the effect of a decreasing 
/gas towards lower cluster masses is amplified and the trend of 
an increasing / star towards lower masses is weakened (at 7?20o)- 
Figure[8] shows that these trends are still present at smaller radii 
and become progressively stronger towards the center. Inter- 
estingly, these effects partially compensate for the baryon frac- 
tion, fh, which attains a much weaker mass dependence. Com- 
paring the radial profiles of these mass fraction for all clus- 
ters (left panel of Figure [8]) and a sub-sample of the third most 
relaxed clusters (right panel) we virtually do not see any dif- 




Figure 8. Radial median profiles of the mass fractions in gas (/gas, solid), stars (/star, dash-dotted), and baryons (/),, dashed). Each profile is normalized by the 
universal baryon fraction and different colors correspond to different cluster mass ranges (indicated in the legends). The profiles in the left panel include all clusters 
and the right panel shows the median of a sub-sample of the third most relaxed clusters, selected by the lowest third in total kinetic-to-thermal energy ratio, K/U. 
In contrast to our non-equilibrium measures such as clumping or kinetic pressure support, the mass fractions in gas, stars, an baryons are almost independent on the 
clusters' dynamical state. 



ferences. In contrast to our non-equilibrium measures such as 
clumping or kinetic pressure support, the mass fractions in gas, 
stars, and baryons are almost independent on the clusters' dy- 
namical state. The only exception are the central values for 
/star in massive clusters which are higher by 5% in our relaxed- 
cluster sample in comparison to all clusters. 

However, there are two noticeable differences of the mass 
fractions in our AGN feedback model in comparison to those 
in our radiative cooling model (Figure|7]). (1) The baryon frac- 
tions do not sum up to the value of our shock heating-only 
model. There is a deficit of gas mass at Z?200 that decreases 
from 10% at M 200 ^ 10 14 M Q down to 5% at 7 x 10 14 M Q . 
Only for the largest clusters (M200 — 10 15 M Q ) does the baryon 
fraction of the AGN feedback model correspond to th at o f the 
shock heating-only model. As explained in in Section l4~Tl this 
is due to the powerful outflows of gas that deposit considerable 
amounts of gas beyond R200 so that the cosmic mean is not even 
reached by 4/?2oo- The only exception to this are the largest 
clusters whose gravitational potential provide a strong enough 
pull on the gas that was successful in having re-accreted the 
gas. (2) In our AGN feedback model at ^200, /star remains al- 
most independent of cluster mass (in contrast to the values at 
^500)- 

This behavior is very interesting with respect to its implica- 
tions for the redshift evolution of / staI on group scales. In order 
to build up high-mass clusters with low values of /tar// ~ 
0.1 -0.15, their progenitors also had to have these low stellar 
mass fractions in a hierarchically growing universe since the 
amount of baryons locked up in (low-mass) stars remains there 
throughout a Hubble timeQ Since / staI is increasing towards 
lower cluster masses in our radiative cooling model, this im- 

2 There can be return of stellar mass to the gaseous phase by means of stel- 
lar winds which is though t to mostly sustain star for mation in galactic disks 
(Martig & Bournaud 2010; Leitner & Kravtsov 2011). If there is a consider- 
able population of intracluster stars iPuchwein et al. 2010), stellar mass loss 
could be a potential channel of increasing /g as at the expense of /star, Dut tne 
low observed mass fraction of intracluster stars renders it unlikely to be an im- 
portant contribution, especially on the high-mass cluster en d I Lin et al. 2003 ; 
IZibetti et alj2 005: Gonzalez et al. 2005, 2007; Giodini et alJ2009h . 



plies that the groups we observe today, do not reflect the prop- 
erties of groups at high-z (in terms of / star ) that merged to as- 
semble large clusters today. This is demonstrated by the almost 
constant behaviour of / tar as a function of cluster mass at z = 1 
in our radiative cooling model (right-hand panel of Figure [7]). 
The values of / star in groups at z — 1 is almost at the same level 
as /tar for big clusters today. In contrast, the approximate con- 
stancy of /tar with cluster mass in our AGN feedback models 
suggests that groups that are observed today do not need to be 
different from group-scale progenitor systems of clusters today. 

The increasing value of / ta r with time on groups scales also 
suggests an excess star formation at late-times which results in 
too blue colors and too massive stell ar masses at the high -mass 
end in comparison to observations (Siiacki et al. 2007b). We 
note that the distributions of / star and / gas are almost Gaussian, 
resulting in almost identical median and mean. 

4.3. Redshift evolution 

First, we show the redshift evolution of the median profiles of 
the stellar and gas mass fractions in our AGN feedback model 
at fixed cluster mass, 1.4 x 1O 14 M < M 2 oo < 2.6 x 1O 14 M 
(Figure [9j. To address the cluster-to-cluster scatter, we ad- 
ditionally show the redshift evolution of the median and per- 
centiles containing 68.3% probability centered on the median 
(corresponding to 1-er intervals) for small and big clusters sep- 
arated at the characteristic mass of M 2 no = 3.1 x 1O 14 M (Fig- 
ured)]). 

For a given cluster-mass bin and with increasing redshift, we 
are probing clusters with a higher mass accretion rate on aver- 
age, i.e., dynamically younger systems. Back to z = 1, there is 
almost no evolution in these mass fractions, especially for big- 
ger clusters and larger radii. The lack of evolution in the stellar 
fractions matches optical observations of galaxy clusters that 
find an abundance of red elliptical galaxies with in clusters that 
have been passively evolving since z > 1 (e.g., Sta nford et alJ 
119981) . This population of passive galaxies i s the foundation 
for th e red sequence cluster finding algorithm ( Glad ders & Yed 
2000), which is the standard cluster finding method in the op- 
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Figure 9. In our AGN feedback model, there is little redshift evolution in the 
median profiles of the gas mass fraction, /g as (< r), (solid) and stellar mass frac- 
tion, / s tai (< r) (dash-dotted), when radii are scaled to ^200- Both quantities are 
normalized by the universal baryon fraction and shown as a function of red- 
shift for fixed mass range (different colors). In contrast, the profiles / ga s(< r) 
that assume hydrostatic equilibrium (dashed) show a larger bias at high redshift 
due to the greater kinetic pressure contribution at earlier times. Those probe 
on average dynamically younger objects that are in the process of assembling. 



gas out (Figure|9l see also iMcCarthy et al.ll201 lh . At z > 1.5, 
we run out of statistics for the mass range under considera- 
tion. However, this interpretation is confirmed by looking at 
a smaller cluster mass range at an even higher redshift which 
shows a further increase in the stellar mass fraction. 

This finding is also consistent with the amount of energy in- 
jected by AGN feedback over the simulation time. Approxi- 
mately one third of total energy is delivered in the cluster for- 
mation phases at z > 2 (analogous to high-z quasar-like feed- 
back), another third in the redshift range 1 < z < 2, and the 
final third below z = 1 (analogous to jet/bubble like feedback). 
These fractions depend moderately on the the numerical reso- 
lution; increasing the resolution enables to resolve the growth 
of smaller halos at earlier times and caus es a higher fraction of 
energy injection at higher redshifts (see Batta gha et al.l l2010. 
for a discussion). 

Despite the almost constant evolution of the gas fractions for 
r > R200, there is positive evolution of the biased gas mass frac- 
tion /g as ,HSE at these radii due to the evolution of the kinetic-to- 
thermal pressure ratio when radii are scaled to 1?200 (Figure[9]). 
Most of this apparent evolution can be absorbed by redefining 
R200 as the radius at which the mean interior density equals 
A time s the mean density of the universe (see Figure 21 in 
BBPS1). However, this rescaling would then imply a redshift 
evolution of / gas with smaller values of / gas at higher redshifts. 
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Figure 10. Redshift evolution of the gas mass fractions, fgas(< r), that are 
computed within different radii (Rjoo> ^200. 2i?2oo. an d 4R200) f° r our AGN 
feedback model. We contrast median values for groups (1 X 10 14 M Q < 
M200 < 3.1 X 10 14 M o , blue) and clusters (3.1 X 10 14 M Q < M 200 < 1 x 
1O I5 M0, blue) and find very little redshift evolution for / g as(< r), particu- 
larly for large clusters. While / gas increases with increasing radius, the baryon 
mass fraction does not reach the cosmic mean for r < 4i?2oo since the stellar 
mass fraction is / s tai (< 4R200) — 7%. 

tical for z < 1 and has been extended to higher redshifts by 
the IRAC Shallow Cluster Survey (ISCS. lBrodwin et al.ll2006h 
and The Spitze r Adaptation of the Red-sequence Cluster Sur- 
vey (SpARCS. iMuzzin et all 120091) . We demonstrate in Fig- 
ure[lO]that there is weak evolution (at the level of 1-cr) for our 
small clusters at /?20o which disappears for larger radii or larger 
clusters. In particular for clusters withM2oo > 3.1 x 1O 14 M , 
we see no evolution of the median / gas already for radii as small 
as R 5m . 

At z = 1.5, both / s tar and / gas are increased towards the center 
indicating that those redshifts probe the tail end of the epoch of 
our AGN feedback which is most efficient at z > 2 in blowing 



4.4. Comparison to data 

In Figure [TT] we compare our AGN feedback model with 
data. On the left-hand panels, we compare radial / gas profiles 
of clusters that are sorted in two bins according to their spec- 
troscopic temperatures, 2.5 < r spec /keV < 5.0 and 5.0keV < 
r spra> We use the r spec - M500 relation by Vik hlinin et alJ 
(2009a) to translate our M500 values to spectroscopic temper- 
atures. While the data in the two upper panels em ploy Chan- 
dra data dVikhlinin et alj|2006t lAllen et al.ll2008l). the bottom 
panel uses data from XMM-Newton dPratt et all 1201 Oh . There 
is generally good agreement of our model clusters populating 
the high-temperature bin, especially with the Chandra data. In 
contrast, our simulated small-mass clusters tend to have on av- 
erage slightly too high / gas values. 

In the right-hand panels of FigureQT| we compare the cluster 
mass dependence of enclosed gas mass fractions (top), stellar 
mass fractions (middle), and baryon mass fractions (bottom). 
Note that we only use those compilations of ft, versus mass 
where the gas and stellar masses have been obtained from the 
same clusters, respectively. Overall, our AGN feedback model 
successfully reproduces observed values and trends, albeit with 
slightly higher values of / gas (and to a lesser extent ft,) at the 
mass scales of groups. While the scatter of our simulated distri- 
bution is noticeably smaller, we caution that we did not include 
measurement biases that would substantially add to the simu- 
lated scatter as demonstrated in Section[3] 

There could be several explanations for an overshooting of 
our simulated / gas values at group scales. Those could in- 
clude numerical resolution, shortcomings of the AGN model- 
ing, or missing physics in the simulations. Increasing resolu- 
tion would enable us to resolve the formation of smaller ob- 
jects earlier on which may allow a more vigorous blow-out of 
the gas already in halos in the mass range 10 12 < M2OO/M0 < 
10 13 . Those are just above the knee of the galaxy luminosity 
function and hence have a decreasing efficiency of star for- 
mation per unit halos mass in more mass ive halos which is 
usually attributed to AGN feedback (e.g.. iBower et aUHOlU: 
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Data comparison of radial f gas profiles: 



Data comparison of enclosed mass fractions: 
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Figure 11. Comparison of our AGN feedback model with data. Left panels : we compare radial fm s profiles of two cluster temperature bins that are separated by 
spectroscopi c temperature r spe c = 5 keV (using the r S pec ~ M500 relat ion by [Vikhlinin et al. 1 2009a) for translating our M500 values). The data samples are taken 
from Vikhlinin et al. 12006), Allen et al. 12008), and Pratt et al. (2010) (top to bottom). Right panels: we compare the cluster mass dependence of enclosed gas mass 
fractions (top), stell ar mass fractions (middl e ), and consisten t treatments of the baryon mass frac ti ons, i.e., /i, obtained f rom the same obj ects, respectively (b ottom). 
We us e /gas data by Sanderson et al. (2003); Vikhlinin et al. (2006, 2009a); Arnaud et al. ( 2007); iGiodini et all <2"0Tj9T) ; ISun etal] f2"00"9T) and / star data bv lLinetalJ 
( 2003); Gonzalez et al. 12007); Giodini et al. ( 2009); Leauthaud et al. ( 201 1). Overall, our AGN feedback model successfully reproduces observed values and trends, 
albeit with slightly higher values of /g as at the mass scales of groups. While the scatter of our simulated distribution is noticeably smaller, we caution that we did not 
include measurement biases that would substantially add to the simulated scatter. 
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iPuchwein & Springelll2012l) . 

As laid out in Section[2] we use an extremely low-resolution 
incarnation of AGN feedback that enables us to simulate large 
volumes and cluster samples, which are needed to provide 
precise predictions for the next generation X-ray and SZ sur- 
veys. Our AGN feedback prescription couples the thermal 
energy feedback to the global SFR of the central cluster re- 
gion. While this model is able to arrest overcooling in mas- 
sive clusters, it may be too simplified for groups. The observed 
(low-frequency) radio emission that spatially coincides with X- 
ray cavities suggests the presence of non-thermal components 
composed of cosmic rays and magnetic fields in these lobes. 
Modellin g AGN feedback i nto cosmic ray rather than thermal 
energy (Siia cki et alj|2008h or exploring kinetic j et feedback 
that thermalizes at larger radii (Du bois et alj |2010) may be vi- 
able alternatives to our simulated scenario, but would also re- 
quire better spatial resolution. 

Finally, our model may be missing important physics that 
could additionally add ent ropy to the centers of groups, such 
as the recently sugges ted (Brode rick et alj|2012t IChang et alj 
l2012tlPfrommer et al.ll2Q12h IPuchwein et alj2012l) blazar heat- 
ing model that proposes the bolometric luminosity of TeV 
blazars as an additional heating source of the inter-galactic 
medium. Phenomenologically, highly biased regions that first 
turn-around to form groups and clusters should be heated first. 
This would imply an evolving pre-heated entropy floor and 
cause an additional decrease of / gas in late-forming groups that 
are observed on a sufficiently short time scale after formation 
so that the central group medium had not have t ime to enter the 
strongly cooling regime (Pfrommer et al. 20T21). 

4.5. Comparison to previous work 

Given the importance of cluster X-ray observations for cos- 
mological parameter estimation, there has been a comparably 
large body of work concerning gas fractions in cosmological 
cluster simulations. Early work on ft, in non-radiative simu- 
lations showed a strong systematic discrepancy between dif- 
ferent numerical methods. This is surprising since all codes 
should solve the same equations for the gas and collision- 
less DM physics without any ambiguity in formulating a sub- 
grid description of star formation and associated feedback pro- 
cesses. SPH simulations of ten massive clusters in a ACDM 
cosmology derived values of ft, ~ 0.87 and 0.83 within R v i r 
and 0.5R vlr w ^500, respectivel y, and show onl y a very weak 
evolution of ft, from z = 1 to dEke et aL| [l998). A number of 
further studies that all e mploy the energy formulation of SPH 
confirm these findings dFrenket all 1 19991: iBialek et alj|2001b 
Mua nwong et alj 120021) . In contrast, cluster studies that em- 
ploy the adaptive mesh refinement (AMR) technique for solv- 
ing the hydrodynamic equations find higher values of ft, « 
0.97 ± 0.03 and 0.94 ± 0.03 with in R yn and 0.5fi vir « #500, 
respectively dKravtsov et al J 120051) . Cluster simulations with 
the entropy-conserving SPH code GADGET yield values of 
ft, « 0.92-0.93, higher by about 5 percentage points than the 
energy-conserving form ulation of SPH but still smaller t han the 
AMR results (see lKav et"aTl 2004; Kr avtsov et alj|2005l which 
is also in line with our findings). All the different methods find 
no significant redshift evolution of the baryon fraction in non- 
radiative simulations. 

Simulations with radiative cooling and star formation nec- 
essarily yield larger discrepancies in the resulting ft, values 
owing to the larger uncertainties of implementing subgrid star 
formation and feedback processes. Common to all of these 
approaches is that too many stars condense out of the ICM, 



which results in the "overcooling problem", i.e., / gas values 
that are too small and / star values that are too high in compar- 
ison to observations. SPH simulations (employing the energy 
formulation) with radiative cooling or feedback yield values of 
ft, ~ 0.85 - 0.9 for the cluster mas s range considered here (see 
Figure 3 of Muanwong et al. 2002). For smaller galaxy clusters 
(within a factor of two of mass M 2 qo = lO'yMp ), entropy- 
conserving SP H simulations by iKay et alj d2004[) and AMR 
simulations by iRravtsov et alj d2005l) find similar average ra- 
dial profiles of / gas . However, those results differ for larger 
cluster s and the AMR re sults obtain substantially higher / gas 
values. lEttori et al] d2006l) presented an analysis for a sample of 
galaxy clusters simulated with the entropy-conserving GAD- 
GET SPH code that followed radiative cooling and star forma- 
tion (essentially similar to our radiative cooling model). They 
find total baryon fractio ns of fb(<Ryii) ~ 0.93 -0.95, similar to 
our values and those by IRravtsov et alj d2005l) on clu ster-mass 
scales. The gas fractions found inlEttori et alj Q2006) and more 
recently in Sembolini et al. (2012) of 0.75-0.8 agree with our 
values (for the radiative cooling model), but are higher than the 
corresponding values bv lKravtsov et alj d2005l) because of the 
different implementations of cooling and star formation. 

More recently, baryon fractions have b een studied in AGN 
feedback models. Puc hwein et alj (|2008) presented mock X- 
ray observations of a sample of zoomed cluster simulations 
which account for radiatively cooling gas of primordial com- 
position, star formation, and thermal energy feedback by black 
holes that grow through Bondi-Hoyle accretion. These simu- 
lations produce gas mass fractions that match the observations 
and are also in agreement with our AGN feedback results for 
larger clusters (M500 >2x 1O 14 M ). However, our smaller 
systems appear to have slightly larger values of / gas due to the 
coarser resolution of our cosmological simulations in compar- 
ison to those zoomed simulations that are able to resolve AGN 
feedback at earlier times in galaxies above L*. In contrast, 
our stellar mass fractions on cluster scales match the observed 
values, but are smaller by about a fact or of three in compari- 
son to the AGN feedback simulation by lPuchwein et alJ (12010) 
(see their Figure 5). Imp roving the simulated physics of these 
AGN feedback models, Fabian et alJ (120101) additionally ac- 
count for metallicity-dependent radiative cooling and follow 
the detailed chemical evolution of the IGM. On cluster-mass 
scales, t hey find slightly small er gas fractions than ours and 
those by Puchwein et al. (2008). Additi onally, they have sim - 
ilarly large stellar mass fractions as IPuchwein et al.1 d2010l) . 
which appear to be in conflict with the data. Common to all that 
body of work (including ours) is that it estimates the black hole 
accretion rate and hence the amount of feedback energy di- 
rectly from simulated hydrodynamical quantities (i.e., entropy 
or density) by means of a subgrid accretion model, which en- 
sures asubstantial predictive power of this approach. In con- 
trast, I Young et al.1 (12011) der ive their feedback energy from a 
semi-analytical model of galaxy formation that they couple to 
hydrodynamical simulations. By construction, their approach 
suffers less from resolution effects in comparison to ours and 
produces a strong cluster-mass dependence of / gas (< Ra) for 
various characteristic cluster radii, in agreement with the ob- 
servations, particularly on group scales. However, their mod- 
els appear to have problems in reproducing the larger / gas val- 
ues of some o f the observed / gas profiles in large clus ters with 
r spec > 5keV dVikhlinin et alJl2006tlAlie"n et alBool . 

Summarizing, it appears that there are still systematic differ- 
ences between Lagrangian SPH and Eulerian AMR codes in 
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computing baryon fraction in non-radiative simulations. Those 
however are reduced to the level of about five percentage points 
when AMR results are compared to those obtained with the 
entropy-conserving GADGET SPH code. Since this is the 
baseline model out of which radiative cooling starts to con- 
dense cool gas and to form stars, more work is needed to 
fully understand the reason of the discrepancy. Purely radiative 
models, which only account for star-formation feedback, suf- 
fer from the over-cooling problem and produce stellar and gas 
mass fractions in conflict with observations. Instead, a growing 
body of AGN feedback models demonstrates that the overcool- 
ing problem can be weakened or even circumvented, depending 
on the detailed implementation of cooling and feedback. More- 
over, since that type of physics causes larger changes in ft, and 
/gas in comparison to the systematic differences owing to the 
numerical method, this holds the promise that the simulations 
are starting to produce robust and predictive results (for a cer- 
tain implementation of AGN feedback) so that a comparison to 
observations is a meaningful and well-posed task. 

5. CONCLUSIONS 

In this paper, we performed a study of gas and stellar mass 
fractions in clusters with the goal to address (1) the bias of 
/gas measurements in recent Suzaku X-ray measurements of 
the ICM beyond the virial radius and (2) provide a theoretical 
framework for using gas masses as a tightly correlated and un- 
biased proxy for the halo masses. This property lies at the heart 
of using X-ray-inferred gas masses to derive precision cosmo- 
logical parameters and is of critical importance for SZ surveys 
as / eas sets the normalization of the Y -M scaling relation and 
the SZ power spectrum scales with the square of / gas . 

5.1. Biases in X-ray inferred gas mass fractions 

We identified two main sources of bias for X-ray inferred 
values of / gas in addition to sample variance of clusters and 
angular orientation variance within clusters; all of which can 
be addressed with cosmological simulations of a large statis- 
tical sample of clusters. (1) Adopting hydrostatic equilibrium 
masses, i.e., neglecting the non-thermal pressure, biases the in- 
ferred masses low by 20-25% on average (or 15-20% if we 
considered only a sample of the third most relaxed clusters). 
Hence / gas is biased high by the same factor that increases dra- 
matically for even larger radii due to the larger relative con- 
tribution of kinetic pressure. (2) The presence of gas density 
clumping biases the X-ray surface brightness and hence the in- 
ferred density and gas masses high. Hence, the resulting values 
for / gas are also biased high at a level of ~ 10-20% within /?2oo 
(slightly mass-dependent and at the low-end for relaxed clus- 
ters). For the most massive clusters with M200 >3x 10 14 M Q , 
it increases even more steeply outside R2Q0 in comparison with 
the bias by HSE masses. 

While these two processes represent an upward measurement 
bias of / gas in comparison to its cosmic value at Z?200> we also 
consider the sample variance of the true / gas as well as its bi- 
ased measurement across clusters and within clusters for differ- 
ent angular directions. Interestingly, the cluster-to-cluster scat- 
ter of the true / gas decreases for larger radii to values around 5% 
at 7?2oo- In contrast, the cluster sample variance of the biased 
measurement of / gas due to the assumption of HSE masses and 
neglecting density clumping remains fairly constant as a func- 
tion of radius at a level of oy„ a ,// ga s — 0.1-0.2 (depending on 
the dynamical state and mass of clusters) within R200 and even 
starts to increase for larger radii due to the strong clumping 



term. 

These different measurement biases are all manifestations of 
two underlying physical processes which both lead to signifi- 
cant angular variance of the (gas and total) mass distributions. 
First, there is a strong internal baryon-to-dark-matter density 
bias, a consequence of collisional-to-collisionless physics as a 
result of the differing efficiencies of gas thermalization and the 
relaxation processes of the DM component. Second, hierarchi- 
cal structure formation implies that cluster are the latest object 
that had have time to form until now which manifests itself in 
a considerable anisotropy of the mass distributions, in particu- 
lar around the virial radius where the filamentary cosmic web 
connects to the cluster interiors. The accretion shock of gas 
that accretes along filaments forms further inwards, leaving a 
larger level of clumping and kinetic pressure contribution al- 
ready around /?2oo- Simultaneously, this causes a considerable 
measure of ellipticity (as measured by the eigen values o f the 
moment-of-inertia tensor) throughout the clusters (BBPS1). 

Equivalently, computing the mass profiles in 48 cones, 
whose footprints partition the sphere, we find anisotropic gas 
and total mass distributions that imply an angular variance of 
/ g as at the level of 30% (almost independent of radius and clus- 
ter mass). Since the angular variance as well as the outliers 
of the gas and DM mass distribution in angular cones is larger 
than the angular variance of / gas , this implies spatial correla- 
tions of the DM and gas distribution through self-bound sub- 
structures that are still holding on to (some) of their gas, es- 
pecially outside /?200- We also demonstrated that in the most 
extreme cases, the cone profiles of / gas can be biased high by a 
factor of two in massive clusters (M200 ~ 10 15 M Q ). While this 
factor will be lower in projection, hydrostatic mass bias as well 
as density clumping may (partially) compensate for these pro- 
jection effects. If the anisotropic mass distribution is indeed re- 
sponsible for some of the high / gas values inferred from Suzaku 
measurements along individual radial arms, it predicts a signif- 
icantly lower level of / gas along other directions. This predic- 
tion can be tested by increasing the number of radial arms in a 
given cluster. 

We caution that the magnitude of the angular variance of the 
mass distributions is not independent of the cluster variance, 
and hence is not additive. Instead, both variances are comple- 
mentary ways of characterizing the non-equilibrium processes 
in clusters (anisotropy, kinetic pressure contribution, clump- 
ing) that become increasingly complicated at the virial radius 
and beyond (BBPS4). We also emphasize that the distributions 
of mass profiles in angular cones is positively skewed (there 
are more extreme positive outliers in M gas and M tot cone pro- 
files) which is expected from the presence of substructures. In 
particular, the outliers with high values of the mass profiles in 
individual cones likely correspond to f ilament s that channel the 
accreting substructures to the cluster (BBPSJ). 

Alternatively the increased level of X-ray surface brightness 
towards the virial region may be due to instrumental or ob- 
servational effects, such as of unresolved X-ray point sources, 
point-spread function leakage from masked point sources, or 
stray light from the bright inner cluster core. 

5.2. Gas mass fraction in cosmological simulations 

We find little redshift dependence of the radial profile of / gas 
in our AGN feedback model, if radii are scaled to i?2oo- In par- 
ticular, this is the case for big clusters withM2oo ^3x 10 14 M Q 
or/and for larger containment radii r > /?5oo- The redshift in- 
dependence of / g as depends critically on the ability to scale the 
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radial variable with Z?200 as the / gas profile increases steadily 
with radius. If the redshift-dependent systematics of inferring 
R200 can be controlled, this provides a good theoretical basis 
of using the gas mass fraction to do precision cosmology. The 
physical reason behind this argument is the deep gravitational 
potential of big clusters (M200 >3x 1O 14 M ). This makes 
it impossible for AGNs to expel gas today, in particular as 
the bubble enthalpy is poorly coupled to the thermal energy 
of the surrounding ICM which results in little heating and en- 
tropy generation of the central ICM on a buoyancy time scale 
(Pfromm er et alj2012l) . Moreover, gas that was pushed beyond 
the virial radius by AGN feedback early-on gets pulled in dur- 
ing assembly of these big systems. 

Cluster mergers may transfer energy from the DM compo- 
nent to the gas, potentially reducing / gas at the virial radius 
from its cosmic value. While our entropy-conserving SPH sim- 
ulations suggest a value of / gas ~ 0.93 when scaled with the 
cosmic value, AMR simulations find somewhat larger values 
of / gas ~ 0.97 (iKravtsov et all 120051) . Radiative physics fur- 
thermore reduces / gas by locking baryons up into stars, causing 
a steeper slope of / gas towards the center. AGN feedback (or 
its quasar-like appearance) is most effective early-on during 
the assem bly time of groups at z > 2 and expels gas far be- 
yond /?2oo (McCar thy et al.1 12.01 1). This causes a much steeper 
dependence of / gas on cluster mass. At the same time, AGN 
feedback (in the incarnation that we modeled) is able to arrest 
overcooling and to drop the stellar mass by a factor of twoQ 
As a result, our AGN feedback model keeps / star almost con- 
stant as a function of redshift (for z < 1) and mass, in particu- 
lar for larger systems. This together implies a significantly in- 
creasing baryon fraction, /baiy, with cluster mass from 0.85 (at 
M 2 oo ^ 1O 14 M ) to 0.94 (atM 200 =i 1O 15 M ). Those scalings 
of / gas and /b ary need to be taken into account (e.g., through a 
marginalization process) when using X-ray and SZ surveys of 
cluster to constrain cosmological parameters. We emphasize 
that AGN feedback produces large effects on various thermo- 
dynamic quantities, including significant changes of stellar and 
gas mass fractions. While the qualitative results of our AGN 
feedback model are in line with approaches by other groups 
and produces robust predictions of the interesting parameter 
space preferred by observations, the detailed scalings and pro- 
files may, however, depend on the particular model realization 
and could vary for different implementations. 

We thank Volker Springel for enlightening discussions and 
careful reading the manuscript. Research in Canada is sup- 
ported by NSERC and CIFAR. Simulations were run on 
SCINET and CITA's Sunnyvale high-performance computing 
clusters. SCINET is funded and supported by CFI, NSERC, 
Ontario, ORF-RE and UofT deans. CP gratefully acknowl- 
edges financial support of the Klaus Tschira Foundation. We 
also thank KITP for their hospitality during the 2011 galaxy 
cluster workshop. KITP is supported by National Science 
Foundation under Grant No. NSF PHY05-51 164. 

REFERENCES 

Akamatsu, H., Hoshino, A., Ishisaki, Y., Ohashi, T., Sato, K., Takei, Y., & Ota, 

N. 2011, PASJ, 63, 1019 
Allen, S. W., Rapetti, D. A., Schmidt, R. W., Ebeling, H., Morris, R. G., & 

Fabian, A. C. 2008, MNRAS, 383, 879 

3 This suppression factor would be larger if we had run at higher resolution 
since models without AGN feedback do not produce numerically converged 
results (at the resolution that we could afford running our cosmological boxes). 



Allen, S. W., Schmidt, R. W., Ebeling, H., Fabian, A. C, & van Speybroeck, 

L. 2004, MNRAS, 353, 457 
Allen, S. W., Schmidt, R. W., & Fabian, A. C. 2002, MNRAS, 334, LI 1 
Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2007, A&A, 474, L37 
Battaglia, N., Bond, J. R., Pfrommer, C, & Sievers, J. L. 2011a, 

arXiv: 1109.3709 
— . 2011b, arXiv:1109.3711 
— . 2012a, in prep. 
— . 2012b, in prep. 

Battaglia, N., Bond, J. R., Pfrommer, C, Sievers, J. L., & Sijacki, D. 2010, 
ApJ, 725, 91 

Battye, R. A., & Weller, J. 2003, Phys. Rev. D, 68, 083506 

Bautz, M. W. et al. 2009, PASJ, 61,1117 

Bialek, J. J., Evrard, A. E., & Mohr, J. J. 2001, ApJ, 555, 597 

Booth, C. M., & Schaye, J. 2009, MNRAS, 398, 53 

Borgani, S., & Kravtsov, A. 2009, arXiv:0906.4370 

Bower, R. G, Benson, A. J., & Crain, R. A. 2012, MNRAS, 422, 2816 

Broderick, A. E., Chang, P., & Pfrommer, C. 2012, ApJ, 752, 22 

Brodwin, M. et al. 2006, ApJ, 651, 791 

Chang, P., Broderick, A. E., & Pfrommer, C. 2012, ApJ, 752, 23 
Churazov, E. et al. 2012, MNRAS, 421, 1123 

Dubois, Y, Devriendt, J., Slyz, A., & Teyssier, R. 2010, MNRAS, 409, 985 

Eckert, D. et al. 2011, arXiv: 1 1 1 1 .0020 

Eke, V. R., Navarro, J. E, & Frenk, C. S. 1998, ApJ, 503, 569 

EnBlin, T. A., Pfrommer, C, Springel, V., & Jubelgas, M. 2007, A&A, 473, 41 

Ettori, S., Dolag, K., Borgani, S., & Murante, G. 2006, MNRAS, 365, 1021 

Ettori, S., Morandi, A., Tozzi, P., Balestra, I., Borgani, S., Rosati, P., Lovisari, 

L., & Terenziani, F. 2009, A&A, 501, 61 
Evrard, A. E. 1990, ApJ, 363, 349 

Fabjan, D., Borgani, S., Tomatore, L., Saro, A., Murante, G, & Dolag, K. 

2010, MNRAS, 401, 1670 

Frenk, C. S. et al. 1999, ApJ, 525, 554 

Galli, S., Bartlett, J. G, & Melchiorri, A. 2012, arXiv:1203.5538 

Gao, L., Navarro, J. E, Frenk, C. S., Jenkins, A., Springel, V., & White, 

S. D. M. 2012, arXiv:1201. 1940 
George, M. R., Fabian, A. C, Sanders, J. S., Young, A. J., & Russell, H. R. 

2009, MNRAS, 395, 657 
Giodini, S. et al. 2009, ApJ, 703, 982 
Gladders, M. D., & Yee, H. K. C. 2000, AJ, 120, 2148 
Gonzalez, A. H., Zabludoff, A. I., & Zaritsky, D. 2005, ApJ, 618, 195 
Gonzalez, A. H., Zaritsky, D., & Zabludoff, A. I. 2007, ApJ, 666, 147 
Gorski, K. M., Hivon, E., Banday, A. J., Wandelt, B. D., Hansen, F. K., 

Reinecke, M., & Bartelmann, M. 2005, ApJ, 622, 759 
Hoshino, A. et al. 2010, PASJ, 62, 371 
Hu, W. 2003, Phys. Rev. D, 67, 081304 

Humphrey, P. J., Buote, D. A., Brighenti, E, Flohic, H. M. L. G, Gastaldello, 

F, & Mathews, W. G. 2012, ApJ, 748, 1 1 
Jubelgas, M., Springel, V., EnBlin, T., & Pfrommer, C. 2008, A&A, 481, 33 
Kawaharada, M. et al. 2010, ApJ, 714, 423 

Kay, S. T., Thomas, P. A., Jenkins, A., & Pearce, F. R. 2004, MNRAS, 355, 
1091 

Khedekar, S., Majumdar, S., & Das, S. 2010, Phys. Rev. D, 82, 041301 
Komatsu, E. et al. 2011, ApJS, 192, 18 

Kravtsov, A. V, Nagai, D., & Vikhlinin, A. A. 2005, ApJ, 625, 588 

Kravtsov, A. V, Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128 

LaRoque, S. J., Bonamente, M., Carlstrom, J. E., Joy, M. K, Nagai, D., Reese, 

E. D., & Dawson, K. S. 2006, ApJ, 652, 917 
Lau, E. T, Kravtsov, A. V, & Nagai, D. 2009, ApJ, 705, 1 129 
Leauthaud, A., Tinker, J., Behroozi, P. S., Busha, M. T, & Wechsler, R. H. 

2011, ApJ, 738,45 

Leitner, S. N., & Kravtsov, A. V. 2011, ApJ, 734, 48 

Lin, Y.-T, Mohr, J. J., & Stanford, S. A. 2003, ApJ, 591, 749 

Majumdar, S., & Mohr, J. J. 2003, ApJ, 585, 603 

— . 2004, ApJ, 613, 41 

Mantz, A., Allen, S. W., Ebeling, H, Rapetti, D., & Drlica-Wagner, A. 2010a, 

MNRAS, 406, 1773 
Mantz, A., Allen, S. W., & Rapetti, D. 2010b, MNRAS, 406, 1805 
Mantz, A., Allen, S. W., Rapetti, D., & Ebeling, H. 2010c, MNRAS, 406, 1759 
Martig, M., & Bournaud, F. 2010, ApJ, 714, L275 
McCarthy, I. G. et al. 2007, MNRAS, 376, 497 

McCarthy, I. G, Schaye, J., Bower, R. G, Ponman, T. J., Booth, C. M., 

Vecchia, C. D., & Springel, V. 2011, MNRAS, 412, 1965 
McCarthy, I. G. et al. 2010, MNRAS, 406, 822 

Miller, E. D., Bautz, M., George, J., Mushotzky, R., Davis, D., & Henry, J. P. 

2011, arXiv:1112.0034 
Muanwong, O., Thomas, P. A., Kay, S. T, & Pearce, F. R. 2002, MNRAS, 336, 

527 



GAS AND STELLAR MASS FRACTIONS IN GALAXY CLUSTERS 



17 



Muzzin, A. et al. 2009, ApJ, 698, 1934 
Nagai, D., & Lau, E. T. 2011, ApJ, 731, L10 
Pfrommer, C, Chang, P., & Broderick, A. E. 2012, ApJ, 752, 24 
Pfrommer, C, EnBlin, T. A., Springel, V., Jubelgas, M., & Dolag, K. 2007, 
MNRAS, 378, 385 

Pfrommer, C, Springel, V., EnBlin, T. A., & Jubelgas, M. 2006, MNRAS, 367, 
113 

Pratt, G. W. et al. 2010, A&A, 511, A85 

Puchwein, E., Pfrommer, C, Springel, V., Broderick, A. E., & Chang, P. 2012, 

MNRAS, 423, 149 
Puchwein, E., Sijacki, D., & Springel, V. 2008, ApJ, 687, L53 
Puchwein, E., & Springel, V. 2012, ArXiv e-prints 

Puchwein, E., Springel, V, Sijacki, D., & Dolag, K. 2010, MNRAS, 406, 936 
Rapetti, D., Allen, S. W., Mantz, A., & Ebeling, H. 2009, MNRAS, 400, 699 
— . 2010, MNRAS, 406, 1796 

Rasia, E., Tormen, G., & Moscardini, L. 2004, MNRAS, 351, 237 
Reiprich, T. H. et al. 2009, A&A, 501, 899 

Roncarelli, M., Ettori, S., Dolag, K., Moscardini, L., Borgani, S., & Murante, 

G. 2006, MNRAS, 373, 1339 
Sadat, R. et al. 2005, A&A, 437, 31 
Sanders, J. S., & Fabian, A. C. 2012, MNRAS, 421, 726 
Sanderson, A. J. R., Ponman, T. J., Finoguenov, A., Lloyd-Davies, E. J., & 

Markevitch, M. 2003, MNRAS, 340, 989 
Sato, T. et al. 2012, arXiv:1203.1700 

Sembolini, F, Yepes, G, De Petris, M., Gottloeber, S., Lamagna, L., & Comis, 

B. 2012, ArXiv e-prints 
Shaw, L. D., Nagai, D., Bhattacharya, S., & Lau, E. T. 2010, ApJ, 725, 1452 

APPENDIX 

A. HORIZONTALLY SLICED SPH KERNEL INTEGRAL 

To speed up the mass binning of a distribution of SPH par- 
ticles into a spherical shell, we integrate the mass portion of 
the SPH-kernel within that shell. We make the approximation 
that the smoothing length h is much smaller than the radius of 
the shell in question, r s h e ii, so the shell can be thought of as two 
plane-parallel slices. We obtain a 9 th order polynomial fit to the 
integrated SPH mass fraction below some normalized distance 
x, which is constrained to be exact at x = {-1 , 0, 1 }: 
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where x = (r-ro)/h for an SPH particle at position xq. The 
fit is everywhere accurate to better than 2 x 10" 3 of the total 
mass. Clearly, this approximation breaks down near the center 
of the cluster (where the shell radius could become compara- 
ble to the smoothing length). However, the constraint at ±1 
guarantees conservation of mass even in the event the approxi- 
mation r s heii 3> h breaks down. 
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